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Abstract: We study the next-to-minimal supersymmetric standard model (NMSSM) 
with the TeV scale mirage mediation, which is known as a solution for the little hierarchy 
problem in supersymmetry. Our previous study showed that 125 GeV Higgs boson is realized 
with 0(10)% fine-tuning for 1.5 TeV gluino (1 TeV stop) mass. The ^ term could be as 
large as 500 GeV without sacrificing the fine-tuning thanks to a cancellation mechanism. 
The singlet-doublet mixing is suppressed by tan/3. In this paper, we further extend this 
analysis. We argue that approximate scale symmetries play a role behind the suppression of 
the singlet-doublet mixing. They reduce the mixing matrix to a simple form that is useful 
to understand the results of the numerical analysis. We perform a comprehensive analysis 
of the fine-tuning including the singlet sector by introducing a simple formula for the fine- 
tuning measure. This shows that the singlet mass of the least fine-tuning is favored by the 
LEP anomaly for moderate tan /3. We also discuss prospects for the precision measurements 
of the Higgs couplings at LHG and ILG and direct/indirect dark matter searches in the 
model. 
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1 Introduction 

The first run of the LHC has finished successfully. The Higgs boson, the final missing piece 
of the standard model (SM), was discovered by ATLAS [1] and CMS [2] and its mass was 
measured as ruh ~ 125 GeV [3]. Subsequently, its spin, parity [4-7] and couplings to the 
SM bosons [8, 9] were measured and appeared to have the property predicted by the SM 
within the current experimental accuracy. Now any new physics model must predict this 
almost SM-like Higgs boson. On the other hand, no sign of the new particles beyond the 
SM was discovered. This casts a serious doubt on the naturalness of the electroweak (EW) 
symmetry breaking and new physics models based on it. 

In particular, supersymmetry (SUSY) has been one of the most popular physics beyond 
the SM to solve the hierarchy problem. In the context of SUSY, the LHC Run I did not find 
expected superpartners, instead put the lower bound of superpartner masses as about 1.4 
TeV for gluino, 1 TeV for the first/second generation squarks [10-13] and 700 GeV for the 
third generation squarks [14-17] in simplified models. Thus, the little hierarchy between 
the Higgs mass and superpartner masses now becomes manifest. Within the framework of 
the minimal supersymmetric standard model (MSSM), the Z-boson mass, mz, is obtained 
as 

ml - - 2 /^^ , ( 1 - 1 ) 
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where and rn^j^ are the soft SUSY breaking scalar mass squared of the up-sector and 
down-sector Higgs helds and is the supersymmetric mass. For example, in the constrained 
MSSM the radiative corrections on are dominated by the gluino mass M 3 , and obtained 
as ~ Thus, if the gluino and stop masses are of 0(1) TeV, we need hne-tuning 

among and to realize the experimental value of mz- Furthermore, the stop 

mass is required to be of 0(1) TeV or larger to lead to the Higgs mass rrih ~ 125 GeV. If 
we further rely on the naturalness as a guiding principle, we need a mechanism of SUSY 
breaking which reconciles this little hierarchy with the notion of naturalness. 

The mirage mediation is one of the interesting mediation mechanisms of SUSY breaking 
[18-20] (also see [21]). It is a mixture of the modulus mediation [22] and the anomaly 
mediation [23] with a certain ratio. In the mirage mediation, the radiative corrections and 
the anomaly mediation contributions cancel each other at a certain energy scale, where 
SUSY spectrum appears as that of the pure modulus mediation. Such an energy scale is 
called the mirage scale. TeV scale mirage mediation sets this scale at TeV scale and it was 
pointed out that the TeV scale mirage mediation can signihcantly ameliorate the above 
hne-tuning problem of the MSSM [24-26].^ 

Although the TeV scale mirage mediation is an attractive scenario to solve the little 
hierarchy problem, there are two unsatisfactory features in the MSSM. First, a typical mass 
scale of the H-term is the gravitino mass (10 — 100 TeV) and a delicate cancellation between 
different contributions is required to realize the correct EW symmetry breaking [19, 24, 31]. 
Second, 125 GeV Higgs mass is still difficult to achieve by 0(10) % tuning of parameters 
because the A-term for the top Yukawa coupling is hxed by the model, which prevents the 
stop mixing from taking the optimal value known as the maximal mixing. A simple solution 
for these problems is an extension to the next-to-minimal supersymmetric standard model 
(NMSSM) [19]. 

The NMSSM is the minimal extension of the MSSM by adding a singlet superheld 
S [32] (see for review e.g. [33, 34]). Here we impose the Z 3 symmetry, which does not 
allow the /U-term, nHuHd, in the superpotential, where Hu and denote the up and 
down-sector Higgs superhelds. Instead, the term such as XSHuHd is allowed. After the 
scalar component of S develops its vacuum expectation value (VEV), the effective /x-term is 
generated as /x = \{S). In the Z 3 symmetric NMSSM, all of the dimensionful parameters are 
originated from SUSY breaking. Hence, the value of ^ is also obtained by SUSY breaking 
effects. That gives us a solution for the so-called /x-problem [35]. The H-term in the MSSM 
is replaced with the A-term, Ax which has a scale of gaugino mass in the mirage mediation 
and solves the first problem in the MSSM. The NMSSM has an additional Higgs quartic 
coupling, A^, in the scalar potential, which is helpful to increase the tree-level Higgs mass 
compared with one in the MSSM. This ameliorates the second problem if tan (3 is small. 
Also the mixing with the singlet and doublet helps to raise the Higgs mass if the singlet 
is light, although it tightens the constraint on the singlet from the LEP Higgs search [36] 
(For recent studies on the mixing effect in NMSSM, see [37-50]). 

^ Certain relations among the soft SUSY breaking parameters at the unification scale may be useful to 
ameliorate fine-tuning [27-30]. 
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The TeV scale mirage mediation was applied to the NMSSM in Ref. [40]. It was 
found that the fine-tuning is improved in the parameter region to realize nih ~ 125 GeV. 
The effective fj, can be considerably larger than the EW scale without sacrificing the fine- 
tuning due to a cancellation mechanism, although it is assumed to be around the EW 
scale in conventional natural SUSY models [51] The mixing between the singlet and 
doublet Higgs bosons is suppressed. In addition, the Higgs sector as well as the neutralino 
sector has a rich structure compared with one in the MSSM by adding the singlet. That 
leads to interesting aspects as shown in Ref. [40].^ In this paper, we further extend the 
analysis in [40] including new phenomenological studies. We show that approximate scale 
symmetries play an important roll to suppress the singlet-doublet mixing in the NMSSM. 
We discuss the suppression holds as far as k ~ 0 and « rrvj^^. We introduce a 

simple formula for the fine-tuning measures and study the fine-tuning of the EW symmetry 
breaking in detail. We show that the singlet mass of the least fine-tuning is favored by 
the LEP anomaly [36, 40, 54-56] when tan /3 is moderately large. We also study other 
phenomenological aspects such as the precision measurement of the Higgs couplings and 
the dark matter search. 

This paper is organized as follows. In section 2, we give a brief review on the mirage 
mediation, and the TeV scale mirage mediation. We introduce the Z 3 symmetric NMSSM 
in section 3 and apply the TeV scale mirage mediation in section 4. In section 5, we study 
phenomenological aspects such as the mass spectrum, fine-tuning in the EW symmetry 
breaking, Higgs couplings and dark matter. Section 6 is devoted to conclusion and discus¬ 
sion. In Appendix A, we show explicitly initial conditions of soft parameters, which are 
induced through the mirage mediation in the NMSSM. 


2 TeV-scale mirage mediation 

Here we give a brief review on the mirage mediation [18]. The mirage mediation is the mix¬ 
ture of the modulus mediation and the anomaly mediation with a certain ratio, which would 
be determined by the modulus stabilization mechanism and SUSY breaking mechanism. In 
the mirage mediation, the gaugino masses are written by 

MaiMcuT) = Mo + "^bagl (2.1) 

where ga and ba are the gauge couplings and their (3 function coefficients, and 772 - 3/2 denotes 
the gravitino mass. We assume that the initial conditions of our SUSY breaking parameters 
are input at the GUT scale, Mgjjt = 2 x 10^® GeV. The first term, Mq, in the right hand 
side denotes the gaugino mass due to the pure modulus mediation, while the second term 
corresponds to the anomaly mediation contribution. In addition, we can write the soft 
scalar mass rrii of matter fields </>* and the so-called A-terms of (j)'^ cj)^ (j)^ corresponding to the 

^As a different approach in this direction, see [52]. 

®See also Ref. [53]. 
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Yukawa couplings yijk as 


"I3/2 


^ijkiMcuT) = ciijkMo - (7j + + Ik)-^ 


mf^McuT) = CjMg - 7i( 


"^ 3 / 2.2 "^3/2 


Stt^ 


Svr^ 


Mo0i, 
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7i = 2 ^ glC[ 
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e, = ^Y.9lC2 
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2“(<^*) 
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'ijk I ) 


ijfcl j 


ii = Stt" 




( 2 . 2 a) 

( 2 . 2 b) 


(2.3a) 

(2.3b) 


(2.3c) 


Here, 7 * denotes the anomalous dimensions of 0® and denotes the quadratic Casimir 

corresponding to the representation of the matter held </>*. In the right hand side, OjjfcMo 
and qMq denote the A-term and soft scalar masses squared due to the pure modulus 
mediation. These coefficients, Uijk and q, are determined by modulus-dependence of the 
Kahler metric of (/>*, (j)^ and cj)^ as well as Yukawa couplings. Indeed, by using the tree- 
level Kahler metric, the coefficient Cj is explicitly calculated as a fractional number such 
as 0,1,1/2,1/3, etc. [22, 26, 58]. We would have 0 (l/ 87 r^) of corrections on c* due to the 
one-loop corrections on the Kahler metric [59]. Such a correction would be important when 
Ci = 0, but that is model-dependent. Here, we consider the case with 


O'ijk — Cj + Cj + Cfc. 


(2.4) 


Such a relation is often satished for 0(1) of Yukawa couplings in explicit string-derived 
supergravity models [ 22 ].^ 

It is convenient to use the following parameter [19], 


"^ 3/2 

Moln(Mpz/m3/2)’ 


(2.5) 


to represent the ratio of the anomaly mediation to the modulus mediation. Here Mpi is the 
reduced Planck scale. 

The mirage mediation has a very important energy scale, that is, the mirage scale 
dehned by. 


M„ 


Mqut 


( 2 . 6 ) 


The above spectrum of the gaugino masses at Mgut leads to [19], 


Ma(M^ir) = Mo, 


^See also [ 60 ]. 
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at the mirage scale. That is, the anomaly mediation contributions and the radiative cor¬ 
rections cancel each other, and the pure modulus mediation appears at the mirage scale. 
Furthermore, the A-terms and the scalar mass squared also satisfy 

AijkiMmh) = (cj + Cj -I- Cfc)Mo, = CiM^, (2.8) 

if the corresponding Yukawa couplings are small enough or if the following conditions are 
satisfied, 

^ijk Q T Cj T 1, (2.9) 

for non-vanishing Yukawa couplings, yijk [19]. 

When a = 2, the mirage scale is around 1 TeV. Then, the above spectrum 

(2.7) and (2.8) is obtained around the TeV scale. That is the TeV scale mirage mediation 
scenario. In particular, there would appear a large gap between Mq and the scalar mass 
rui with Cj ~ 0. We will apply the TeV scale mirage scenario to the NMSSM in section 4. 

In the TeV scale mirage scenario, the stop mass squared becomes negative at high 
energy [61], while it is positive at low energy below 10® GeV. Thus, the vacuum which 
breaks the EW symmetry at the EW scale might be a local minimum, but instead there 
would be a color and/or charge breaking vacuum with field values larger than 10® GeV. Here, 
we assume the thermal history of the Universe such that field values remain around the 
origin until the EW scale temperature. In addition, we need to confirm that the tunneling 
rate is small enough, i.e. less than the Hubble expansion rate. In Refs. [62], it has been 
shown that such a rate is small enough, as long as the squark/slepton masses squared are 
vanishing or positive around 10^ GeV. This condition is satisfied in our TeV scale mirage 
mediation scenario. 

3 NMSSM 

Here, we briefly review on the NMSSM, in particular its Higgs sector. We extend the 
MSSM by adding a singlet chiral multiplet S and imposing a Z 3 symmetry. Then, the 
superpotential of the Higgs sector is written as 

VFffiggs = -XSHuHd + (3.1) 

Here and hereafter, for S, Hu and Hd we use the convention that the superfield and its 
lowest component are denoted by the same letter. The full superpotential also includes the 
Yukawa coupling terms between the matter fields and the Higgs fields, which are the same 
as those in the MSSM. 

The following soft SUSY breaking terms in the Higgs sector are induced, 

14oft = — XA\SHuH(i + —Ai^S^ -\- h.c. (3.2) 

Then, the scalar potential of the neutral Higgs fields is given as 

U = A2|5|2(|i7®|2 + |F®|2) + _ Ai^XT + Vd 

X-cn\[J\Hu\’^ + + cng\S\’^ — \A\SHuHfi -\- —Af^S^ -|- h.c., (3.3) 
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(3.4) 


VD = \{gl + gl){\H^d?-\Hl\‘^)\ 

where gi and g2 denote the gauge couplings of U(1)y and SU(2). Similarly, there appear 
the soft SUSY breaking terms including squarks and sleptons as well as gaugino masses. 
These are the same as those in the MSSM. 

The minimum of the Higgs potential is obtained by analyzing the stationary conditions 
of the Higgs potential, 

= \^v cos /3sin^ (3) — \kvs^ sin /3 + ^ g^v^ cos /? cos 2/3 
dH^ 4 

+ rnjj^v cos P — XA\vs sin f3 = 0, (3.5a) 

= X^v sin/3(s^ + cos^ (3) — Xkvs^ cos /3 — \g^v^ sin /3 cos 2/3 
Uiiu 4 

+ m^^t;sin/3 — AHAt'scos/3 = 0, (3.5b) 

= X^sv^ + 2 k^s^ — Xkv^s sin 2/3 + m|s — ^XA\v^ sin 2/3 + kA^s'^ = 0, (3.5c) 

uo Z 

where g^ = <3i + Here, we denote VEVs as 

r;2 = |(i/0)|2 + |(i 70 )| 2 ^ tan/3 = |||, s = (5). (3.6) 

Using the above stationary conditions, we obtain the Z boson mass = \g‘^v‘^ as 


ml = 


1 — cos 2/3 
cos 2/3 




1 + cos 2/3 
cos 2/3 






(3.7) 


where g = As. For tan/3 /:!> 1, this equation becomes eq. (1.1). That is, this relation is 
the same as the one in the MSSM. Thus, the natural values of \mH„\ and \g\ would be of 
0(100) GeV. Furthermore, the natural value of {mn^/ tan/3 would be of 0(100) GeV or 
smaller. Alternatively, \g\ and {muj/ tan/3 could be larger than 0(100) GeV when and 
tan^ (3 are canceled each other in the above relation at a certain level. Fven in such a 
case, \mHu\ would be naturally of 0(100) GeV. On the other hand, other sfermion masses 
as well as gaugino masses must be heavy as the recent LHG results suggested. To realize 
such a spectrum, we apply the TeV scale mirage mediation in the next section, where we 
take ch^ = 0 to realize a suppressed value of \m}{J\ compared with Mq. 


4 TeV scale mirage mediation in NMSSM 

Here, we apply the TeV scale mirage mediation scenario to the NMSSM and study its 
phenomenological aspects. 

4.1 Model 

Soft SUSY breaking terms are obtained through the generic formulas (2.1) and (2.2) with 
taking a = 2. For concreteness, we give explicit results of all the soft SUSY breaking terms 
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for the NMSSM in Appendix A. We concentrate on the Higgs sector as well as gauginos 
and stops. 

A concrete model in the mirage mediation is fixed by choosing c*. We consider the 
following values of Cj [40], 


C-Hd 1) 0) 0) Qi Ctfl 2’ (4'^) 

up to one-loop corrections for Hu, S, and left and right-handed (s)top fields, respec¬ 
tively. This is the same assignment as the pattern II in Ref. [26] for the MSSM except for 
cs- Then, the soft parameters due to only modulus mediation contribution are given by 


(^t)modulus — (^A)modulus — ^ 0 , 

(.'^Hd)irLodu\us ~ ^0 1 {jn_ 


(4 


k; J modulus 


= 0 , 




modulus 


= K-J 


modulus 


= -M 




0 > 


(4.2) 


)modulus (^s)modulus Oj 


up to one-loop corrections. The above assignment of Cj (4.1) satisfies the condition, (2.9) 
for the top Yukawa coupling and the coupling A, but not for the coupling k. However, we 
do not consider a large value of k to avoid the blow-up of k and A as will be shown later. 
Thus, we obtain the following values. 


Ai Aa Mo, (4.3) 

^Hd - Ml ml ^ ml ^^Ml 

up to 0{k'^/8tt'^) at the TeV scale. 

Similarly, at the TeV scale we can obtain 

ml^ « 0, « 0, (4.4) 

up to 0{Mq /SttI, and 

4, 0, (4.5) 

up to 0(Mo/87r^). That is, values of lA^P, and are suppressed compared with 
Mq, and their explicit values depend on the one-loop corrections on the Kahler metric. 
Thus, we use A^ as a free parameter, which must be small compared with Mq. In addition, 
we determine the values of and /x (= As) at the EW scale from the stationary 

conditions, (3.5), where we use the experimental value mz = = 91.19 GeV and tan/3 

as a free parameter. 

Through the above procedure, the parameters, and /x, at the EW scale are 
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expressed by tan/3, as follows ^ 


^ = A(5) = 


A\ tan j3 
2 (l — ^ tan/3) 


{l-Vl-4X}, 

of/ 


ml = -2 (f) A." - (f) +^ml\{-A+2l) sin2/S - 2 


2 ^ tan^ /3 - 1 / 

tan2 ^ \ tan2 /3 - 1 




m- 


where, 


X = 


(1- I tan/3) r tan^/3 / 
tan^ /3 I tan^ /3 + 1 y 


2 A^ tan^ (3 — 1 
2 


2 tan^ /3 

For tan/3 max(l,K/A), these parameters are approximated as 


ml 


mu, 




rog ~ —2 


m 




A / I 24a tan (3 


Al, 


+ 2 


A 2 242 


( ^ 

y24Atan/3 j ' ^2^^ 


A ^2 

,2 


m 


mu,. 


H, 


m 




rur 


tan^ (3 A\ tan^ (3 


(4.6a) 

(4.6b) 

(4.6c) 


(4.7) 


(4.8a) 


(4.8b) 


(4.8c) 


These formulas tell that the leading terms in the expansion by tan/3 are O(m^) for nig 
and and the next-to-leading terms are 0{Mq / tan^ (3) because ~ 24 a — Mq. 

If tan/3 = 0(10), the values of g and could be of 0(100) GeV while the other- 

masses of the superpartners are of 0{Mq) = 0(1) TeV. Thus, the hne-tuning problem 
can be ameliorated. Actually, the hrst and the second terms in the last equation cancel 
each other for our choice of c^. The next leading contributions are of 0{m^^/tarA (3) 
or 0{m‘jj^g/ taiA (3Ax). Then, ml is almost determined by alone and insensitive 

to ~ Mo/tan/3. This means that tan/3 ~ 3 is enough to obtain the hne-tuning of 
|cllnm|/(91nm|^^|“^ = ml/2m‘jj^ = 0(100)% for Mq ~ 1 TeV. Therefore g can be as 
heavy as 0(400) GeV without sacrihcing the hne-tuning. This cancellation originates in 
the structure of the doublet mass matrix. 


CM = -{Hd,K)Mjj 



(4.9) 


where 


Ml = 


mjj^ + lA 


-Axg 


Mq2 -a ijA -Mog\ 

-Mog j ■ 


An 777^ -I- ^ 

-Axg m^^ + g 

Galculating the determinant of the mass matrix, it can be easily checked that the modulus 
mediated contribution Mq cancels, det(A4|^) ~ while the trace of the mass matrix is 


®We are only interested in the natural spectrum. Thus we discard another solution of the quadratic 
equation, fj, « Ax tan/3/2. 
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Mq + 2^^. This means that the heavy mode has mass of 0{Mq) and the mass of the 
light mode is suppressed as /r^/Mo ~ tan 13. Therefore a flat direction appears along 
Hu/Hd ~ ~ tan/3. This mechanism was previously observed in [26] in the context 

of the MSSM where Ax is replaced with the S-term. In the mirage mediation, i3-term is 
a remnant of the fine-tuned cancellation between the terms of 0 ( 7713 / 2 )- Thus the relation 
rriHa ~ is subject to uncontrolled corrections. However, in the NMSSM, the relation is 
well controlled up to the leading contribution of the modulus mediation. 

This cancellation is also important for raising the SM-like Higgs boson mass radiatively. 
In the NMSSM, singlino-higgsino loop with A = 0(1) could contribute to the radiative 
correction of the Higgs boson mass in addition to the top quark loop [41, 63]. If the singlet 
scalar is as light as the singlino or higgsino, this correction is canceled by the scalar loops 
because of SUSY. However, the structure of the mass matrix (4.10) tells the coupling \S\'^h'^ 
(~ vanishes due to the SUSY breaking effect. This means the correction could be 

sizable even if the singlet scalar is light as in our case. The renormalization group equation 
in the effective theory roughly gives 


/\rn, ~ - 

^ ~ 167r2 


_4A2!I^ + 8A4-4A^ 

yZ 




- 16A^ 




m 


H 


1 - 


^2 


167r2 


,mt 


-4A^^ + 8A^ 


m 


In 


1 - 

m 

1 


H 


tan2/3 
/ niH 
\^h 


In 


/ ^ 
V '^h 


~ (56GeV)^A2 (A^ - 0.26) (l +1.1 log 


Ho 


rriH 

LlTeVJ 


(4.11) 


where nih ( 777 . 77 ) is the mass of the light (heavy) doublet. In our case, uih ~ ~ Mq. 

This correction amounts ~ 5 GeV for ruh = rnz and A = 0.7. 

Potentially, the mixing between the light doublet and the singlet may invalidate the 
above discussion on the hierarchy of the Higgs boson masses. In models with m'g ~ 0 and 
7C << 1, this mixing is automatically suppressed by \v/Ax/ tan P if we neglect the gauge 
interaction (This suppression is noticed in [37] without a reason.). In this limit, the Higgs 
action has the following approximate scaling symmetry: 

Hu{x) = (4.12a) 

Hd{x) = e‘^^H'aie^x) (4.12b) 

S{x) = S'{e‘^x), (4.12c) 


where (/> is a free scaling parameter. This symmetry is explicitly broken by the Kahler 
potential, 

}Cs = S^S, (4.13) 

the D-term potential and all the kinetic terms. After the Higgs bosons develop the vacuum 
expectation values, the symmetry is spontaneously broken and the Nambu-Goldstone boson 
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appears, which corresponds to the light doublet®. Thus it is inherently mass eigenstate 
without the singlet component. The symmetry breaking Kahler potential generates the 
F-term potential: 

Vf = X^\H^\^\Hd\^ (4.14) 

which does not affect the singlet sector at all. This potential gives additional terms in the 
doublet mass matrix which we neglected before: 


d'^V^ 

dH*dHu 




d‘^V§ 

dH*dHd 




d‘^V§ 

dH*dHu 


\ ^ ZJ* u 

A 11.,Ilf 


d- 


(4.15) 


The contributions to the diagonal elements are eliminated by the minimization conditions 
(3.5a), (3.5b). The CP even part of the mass matrix is given by 


Ml = 


( 


dct)id(l)j 

^2 sin2 /3 


{—A\ + 2X^v‘^) sin (3 cos /3 XA^v cos 2/3 sin /3 


(—^2 + 2X^v‘^) sin /3 cos /3 
cos 2/3 sin /3 


V 

( cos /3 — sin /3 0 
sin /3 cos /3 0 
0 0 1 


^2 cos2 /3 

—XAxv cos 2/3 cos j3 


( 0 
0 


0 


y 0 —XA\v cos 2/3 


0 

—XA\v cos 2/3 
X^v^ 


+ 


/ 0 A2z;2sin2/3 o' 

A2?;2sin2/3 0 0 

0 0 0 


—XAxv cos 2/3 cos /3 
A 2^;2 

^ cos /3 sin (3 0 
— sin /3 cos /3 0 

y 0 0 1 


(4.16) 


where (/>* = x/2{ReHd — {Hd), ReH^ — {Hu), ReS — (S)) and we used the minimization 
condition (3.5c), e.g. = ^Aa sin 2/3 which is not affected by the gauge coupling. The 
last term which corresponds to the off-diagonal contribution in (4.15) mixes the light and 
heavy doublets and slightly changes the diagonalization angle of the CP even mass matrix, 
[3' = 7t/ 2 — a from the direction of the Nambu-Goldstone boson, [3 as 


I 

- tan 2(3' 


- tan 2/3 



25(3 _ _2^ 

sin 2(3 cos 2/3 A\ ’ 


(4.17) 


where 6(3 = (3' — (3. This small rotation generates a non-zero entry in the singlet-light 
doublet mixing terms (1-3, 3-1 elements) in the CP even Higgs mass matrix. The 2-3 (3-2) 
element is —XAxv cos2(3. Thus the rotation yields ~ 2 A^u^/Ha sin 2/3 cos2 2/3 in 1-3 (3-1) 
element. The mass of the light doublet and that of the singlet are order of Xv. Then the 
mixing angle between the light doublet and the singlet is order of Xv/Ax/ tan/3. Similar 

^Similar symmetry can be defined in MSSM, which forbids all the quartic potential. In such a case, the 
non-trivial VEV is not possible unless the mass parameters are tuned so that a massless mode appears in 
that direction. 
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discussion applies for the D-term potential. The diagonalization matrix 5 of is given 

by 


c'f = 


/ cos^ 2/3 + sin^ 2/3 


0 


0 


V 

+0 

( 


= 


t)' 

cos 13 
sin j3 


— sin (3 
cos j3 


Af + (m| - A^v^) sin2 2/3 0 

0 A^v^ sin^ 2/3 


Af-sin^ 


A^cos/3 \ + 0{ — 


—A^ sin2;0 —A^ cos2/3 1 


A\ 


(4.18) 

(4.19) 


at tree level. Note that the mass of the light doublet originates in the symmetry breaking 
terms. Non-vanishing k and are expected to add corrections of order k and m‘gfA‘^v‘^ 
to S. Radiative corrections to the Higgs potential are another source of the singlet-doublet 
mixing because the kinetic terms and the quark-lepton Yukawa couplings also break the 
scaling symmetry. Actually, this is not the end of story. The Higgs action has another 
approximate scale symmetry if k = = 0, in which the transformations of Hu and S 

are interchanged in eq.(4.12). The VEVs of S and spontaneously break this symmetry 
and the Nambu-Goldstone boson appears, which consists of a combination of S and 
corresponding to the singlet-like Higgs boson. In the n = ~ i^d) = 0 limit, the 

real parts of S and Hu are the Nambu-Goldstone bosons (mass eigenstates). After we turn 
on the symmetry breaking terms, the 1-3 and 3-2 elements of the mixing matrix S must 
break the both symmetries because S or Hu is still the Nambu-Goldstone boson if one 
of the symmetries is preserved. In addition, the mixing between S and Hu must vanish 
in the Lagrangian, eq.(3.3) if H^ decouples [38]. Therefore the leading contributions to 
( 5 'I') 3 i, (5'^)23 are given by k{S)/Mo, (Ha)/M o, {Hu){S)/M^, {Hu)m\^/M^, {S)ml/M^ 
and rr/jj^m/g/MQ at tree level. Our argument on the suppression of the singlet-doublet 
mixing generally holds for models with small k and 

The GP-odd part of the mass matrix is diagonalized as 


Ml = 




d'^V 

drjidrjj 




/o 0 0 

0 Al + A^u^ 0 

yO 0 0 


^ cos (3 sin (3 cos 7 sin /3 sin 7 
— sin j3 cos (3 cos 7 cos j3 sin 7 
y 0 — sin 7 cos 7 


V, 


(4.20) 


(4.21) 


in the k = 0 limit. Here r/j = ^/2{Im{Hd), Im{Hu), Im{S)) and tan 7 = Av/A\. Two 
Nambu-Goldstone bosons appear for the broken electroweak and Pecci-Quinn symmetries. 
The former is absorbed by Z boson and the latter acquires mass with non-vanishing k. In 
the following section, we show numerically the spectrum of our model. 
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5 Spectrum and phenomenological aspects 


Here, we study numerically the spectrum and phenomenological aspects of our model. 
Our numerical calculation is performed with NMSSMTools package [64], We choose input 
parameters as the squark and slepton masses, their tri-linear couplings, and parameters in 
the Higgs sector (A, k, A\, and tan/3). We hrst minimize the one-loop effective 

potential of our model at the SUSY scale and calculate the effective /r-term. Then use 
it as an input of the NMSSMTools instead of NMSSMTools allows various levels of 

sophistication in the calculation of the Higgs boson mass. To keep the consistency with the 
input, we choose the default option (BLOCK MODSEL 8 O) which corresponds to the one-loop 
calculation at SUSY scale. The small parameters m|, are calculated so that we can 
obtain the correct EW symmetry breaking. In the following, we only adopt the parameters 
which predict |m||, IrrijjJ ^ Mg/Svr^ within the range of the model ambiguity. We choose 
the small parameter as an input because the number of minimization conditions is not 
enough to predict it. 

5.1 Higgs Mass Spectrum 

We hrst discuss the mass spectrum of the Higgs sector. As described in the previous 
section, our model is in the decoupling region where the heavy doublet cos j3Hu — sin 
orthogonal to the hat direction has a mass much larger than the EW scale. Then this heavy 
doublet approximately consists of the heaviest CP-even, CP-odd and Charged Higgs bosons 
and they have degenerate masses ~ Mq up to the effect of the EW symmetry breaking. 
Remaining light degrees of freedom are the hrst and second lightest CP-even Higgs bosons 
and the lightest CP-odd Higgs boson. 

In hgure 1, we show the mass spectrum of the light CP-even Higgs bosons as a function 
of A and n for tan/3 = 3 and Mq = 1500 GeV. In addition we present and ^ 

to check whether the parameters ht in the range favored by the mirage mediation or not. 
We also plot the coupling squared of the lightest CP-even Higgs boson to the gauge boson 
normalized with its SM value. This measures the amount of the doublet component in the 
lightest CP-even Higgs boson. The red curve in the hgure indicates the boundary where A 
and K blow up at the Planck scale. Inside this curve the model remains perturbative until 
the Planck scale. The gray shaded area around the red curve shows the region where our 
iterative method fails to calculate the minimum of the effective potential. This is due to 
the tachyonic tree-level Higgs mass which we use as the initial condition of the iteration. 
The quantum correction may lift it, however, we do not pursue this possibility because the 
region is already excluded by the LEP Higgs mass limit. The gray shaded area in the small 
K shows the region where the lightest CP-even Higgs becomes tachyonic. The yellow region 
is disfavored due to the false vacuum developed far beyond the EW scale as studied in [65]. 

As we see in the upper left panel, the lightest CP-even Higgs boson mass can reach 
125 GeV for large value of A (~ 0.7) and small value of k (~ 0.1). This region is favored 
by the mirage mediation as shown in the middle left and middle right panels where the 
corresponding small mass parameters satisfy rrig, for the observed Higgs 

boson mass. The lower right panel also indicates hi is almost doublet in this region. This 
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Figure 1. The CP-even Higgs masses and the Higgs mass parameters for tan/3 = 3, Mq = 1500 
GeV, = -100 GeV. 


is because the mixing between the doublet and the singlet in the mass matrix reduces the 
smaller eigenvalue while raises the larger eigenvalue [37-50]. Then the mass of the lightest 
Higgs boson is maximized when the mixing is minimized. The second lightest CP-even 
Higgs is around 200 GeV in the upper right panel and dominantly consists of the singlet. 
The higgsino mass parameter fi is around 500 GeV in the lower left panel which is much 
heavier than the EW scale, however, the cancellation mechanism described in the previous 
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Figure 2. The CP-even Higgs mass and the Higgs mass parameters for tan/3 = 10, Mq = 1500 
GeV, = -100 GeV. 


section protects the model from the fine-tuning which we will see numerically later. It is 
noted that the region is naturally realized in the model where A has a strong dynamics origin 
at the Planck scale and approximate Pecci-Quinn (PQ) symmetry moderately suppresses 

K. 

In figure 2, we plot the same observables for tan/3 = 10. The other input parameters 
are fixed as in figure 1. Now the lightest CP-even Higgs boson mass cannot reach 125 GeV. 
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Figure 3. The lightest CP-odd Higgs mass for 
input parameters are same as figure 1. 





/3 = 3 (left) and tan/3 = 10 (right). The other 


Instead, the second lightest CP-even Higgs boson can be as light as 125 GeV. This region 
is favored by the mirage mediation according to the contours of the small parameters, 
and ms- In this region, the coupling squared to the gauge boson indicates that hi almost 
consists of the singlet and hence /12 is almost doublet. Now the mixing between the doublet 
and the singlet in the mass matrix is non-negligible because the radiative correction is not 
enough to push the doublet to 125 GeV and the mixing is required to raise the larger 
eigenvalue. The F-term contribution to the doublet mass from A is suppressed due to the 
large tan /3. This mixing imposes a rather strong constraint on the model from the LEP 
Higgs boson search via e'’'e“ —?■ Z* —)• Zh. The region around mh^ > 100 GeV still survives 
and is an interesting possibility in light of the observed excess around 98 GeV [36, 40, 54, 55]. 
If we raise the SUSY scale, the radiative correction takes the place of the mixing in realizing 
the 125 GeV Higgs boson mass and the constraint is relaxed. The higgsino mass is around 
Mq/ tan/3 ~ 150 GeV as expected. The LEP chargino mass bound implies tan,0 cannot 
exceed ~ Mo/(100GeV). 

For completeness, we present the contours for the lightest CP-odd Higgs boson mass 
in figure 3. The mass decreases as k is reduced as expected from the approximated PQ 
symmetry. The contours are almost similar to those of the second lightest CP-even Higgs 
boson. 125 GeV CP-even Higgs boson corresponds to ma^ — 200 GeV for tan /3 = 3 and 
ruaj ~ 150 GeV for tan/3 = 10. 


5.2 Fine-tuning in the electroweak symmetry breaking 

Here we numerically estimate the degree of fine-tuning of the EW symmetry breaking in 
our model. Following the standard lore [66-70] we define the fine-tuning measure of a 
observable y against an input parameter x as 


^ gln(//) 

^ ain(x)' 


(5.1) 


To evaluate the EW symmetry breaking, the simplest choice of y’s are the three Higgs 
VEVs Vi = {{Hu), {Hd), {S)). Instead of them, we choose tan/3 and //, extending the 
standard choice, y = m? 2 : the literature [56, 71, 72]. Here we emphasize that in principle 
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all the VEVs must be examined to estimate the degree of fine-tuning once we extend the 
EW sector of the SM model, including usually disregarded tan/3. The maximum of these 
measures, A = max|A|| could be regarded as the measure of the EW fine-tuning in the 
model. As x, we take A, k and the small parameters, ms and at the SUSY scale. 

Note that the large parameters such as mn^, Ax are not free parameters in our model and 
fixed by the ultraviolet physics. 

These measures are easily calculated by the potential minimization conditions eq.(3.5) 
[69]. If we vary the input parameters Xa with Ax^, the Higgs VEVs Vi shift with Avi 
following the potential minimum and they are related as 


E 


d^V 

dipidfpk 


{V2Avk) + ^ 

a 


d^V 

dcpidxa 


Axa = 0. 


(5.2) 


Noticing that the coefficient of the first term is given by the CP even Higgs mass matrix 
defined in eq.(4.16), this leads to 


dvj 
cHn Xa 




d‘^V 


d4)kdlnxa 


(5.3) 


The fine-tuning measures can be constructed from these derivatives as 



2 

V 



dvd 

cHnx 


-|- sin /3 


dvu '\ 
dlnx) ’ 


(5.4a) 


^tan /3 _ 




1 dvu 1 dvd 


Xu 5 In X Vddlnx^ 
+ SxX- 


Id dlnx 


(5.4b) 

(5.4c) 


This expression has clear physical meaning b It is obvious that the measures increase 
if a certain Higgs mass is much lighter than the mass parameters in the potential. In 
NMSSMTools, the Higgs masses are calculated including the quantum corrections, however, 
the fine-tuning measures are calculated using the tree-level formula. This overestimates 
the fine-tuning if any of the tree-level masses vanish although the quantum corrections lift 
them. We calculate the measure using the above formula with the tree-level potential V 
and the quantum corrected mass matrix A4|. 


In figure 4, we show the fine-tuning measures, A^, ^ for tan/3 = 3, Mq = 1500 GeV, 
Ak = —100 GeV. The red dots satisfy the Higgs mass condition, 124 GeV < mh^ < 126 GeV. 
We also plot blue dots representing the region where 124 GeV < mh 2 < 126 GeV. Both of 
them are favored by the mirage mediation within the 1-loop uncertainty. Here we choose 
the width, 125 ± 1 GeV as a guide to the eye. Note that currently the theoretical error of the 
SM-like Higgs mass in the NMSSM could reach 5 GeV [57] and actual bound is somewhat 
weaker than the bands shown in the figure The upper left and upper right panels show the 
case for A and k, respectively. Because the effective /r-term is roughly proportional to k/A as 


^Essentially the same formula for y = m% but different presentation is given in the appendix C of [47]. 
®We thank F. Staub for pointing out this issue. 
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Figure 4. The fine-tuning measures of m% tan /3 = 3, Mq = 1500 GeV, = —100 GeV. A and 

K are scanned over the white and yellow regions in figure 1. 


17- 

































1000 


-< 

oc 

o 


100 


“ 10 
■a 


a. 

oc 

o 

-o 


1 


0.1 


0.01 


124 GeV < m(h,) < 126 GeV 
124 GeV < in(hj < 126 GeV 

\ 

''x_ 

1000 

100 

CiC 

o 

"" 10 

-O 

^ 1 

124 GeV < ra(h,) < 126 GeV - 
124 GeV < m(hj) < 126 GeV 

• 

CiC 


- 

o 



•—• 

- 

Mo = 1500 GeV 

0 1 

Mq = 1500 GeV 

A.= -100 GeV 


A.= -100 GeV 

tan|8= 3 


tan/?= 3 

m^= 172.9 GeV 

0.01 

m^= 172.9 GeV 


0 100 200 0 100 200 


m(h,) 


m(h,) 


1000 


100 , 


124 GeV < m(h,) < 128 GeV : 
124 GeV < m(hj) < 126 GeV 

M„ = 1500 GeV I 
A,= -100 GeV : 
tan|8= 3 
m^= 172.9 GeV 



200 


1000 - 


< 


100 i 


124 GeV < m(h,) < 126 GeV 
124 GeV < m(hj < 126 GeV 


10 


a 


0.1 


0.01 


= 1500 GeV 
A,= -100 GeV 
tan0= 3 
m,= 172.9 GeV 


100 

l(h,) 


200 


1000 


100 


O 10 


0.1 


0.01 


124 GeV < ra(h,) < 126 GeV 
124 GeV < m(hj < 126 GeV 


•A' 


M„ = 1500 GeV 
A,= -100 GeV 
tan/S= 3 
m^= 172.9 GeV 


100 
m(h ) 


200 


Figure 5. The fine-tuning measures of /r for tan/3 = 3, Mq = 1500 GeV, = —100 GeV. 


seen in the figure 2, severe fine-tuning is expected for these parameters. However, ^ ^ 5 

TUp' 771^ 

and Ak ^ < 4 are realized for rrih^ ~ 125 GeV in our model. Also A^ ^ < 10 is satisfied 
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Figure 6. The fine-tuning measures of the EW symmetry breaking for tan/3 = 10, Mq = 1500 
GeV, = —100 GeV. A and k are scanned over the white and yellow regions in figure 2. Only the 
points satisfying 124GeV < < 126GeV are depicted. 


for m/i 2 ~ 125 GeV case unless is not too light. This is not trivial because the khaki 
circles which do not necessarily satisfy the conditions required by the mirage mediation 
reach ^ ~ 100 even with a moderate value of rrih^. We confirmed that the fine-tuning 
measures for these dimensionless parameters do not deteriorate even if we increase Mq to 
5 TeV. The middle left and middle right panels present the results for mjj and m|. It is 

?72^ TTl^ 

found that A 9 ^ 5 and A f 1 for the observed Higgs mass and a moderate value of 

mir mi ~ °° 

Mu o 

m? 

Again, this is not trivial because the khaki circles spread over A f 2 ~ 100. These 

Hu ’ S 

measures for the dimensionful parameters are sensitive to the SUSY scale, Mq. The lower 
left panel shows the case for and we find A^^ < 0.4 for mh^ ~ 125 GeV and A^^ < 10 
for mh 2 ~ 125 GeV with a moderate value of ruhi- This measure appears to be not sensitive 
to the SUSY scale, Mq. 

In figure 5, we show the fine-tuning measures for /r with the same input parameters. 
All of them are smaller than 0(1) for the observed Higgs mass. On the other hand, 
« —A^ holds well in our model because of the relation eq.(4.8a). Then do not 

work as independent measures. Thus the worst measure is A™^, which is below 10 unless 
hi is extremely light (such a region is actually excluded by phenomenological constraints). 
Therefore the model requires at most 10% tuning (A“^ = 0.1) for tan /3 = 3 and Mq = 1500 
GeV. In our model the stop mass is ,ycf~^MQ ~ 1 TeV. Thus we numerically confirmed 
the model exhibits surprisingly low level of fine-tuning compared to the conventional SUSY 
breaking models. 

In figure 6 , we summarize the same fine-tuning measures for the tan /3 = 10 case. The 
other input parameters are same as in figure 4. Here, we impose 124 GeV < nih^ < 126 GeV. 
There is no parameter region satisfying ~ 125 GeV in this case. The fine-tuning 

TnP‘ 

measures for A, k and are below 0(1). The worst tuning (A) is A^f ~ 5 for > 50 
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Figure 7. The fine-tuning measures A™f (left) and A ^2 (right) for Mq = 3 TeV and Mg = 5 
TeV (tan/3 = 20) compared with Mg = 1.5 TeV (tan/3 = 10) case. is fixed at = —100 GeV. 

2 

GeV. While A^f, 2 and A^ 2 are well below the unity with m/i, ~ 100 GeV, they 

77T Q TTl/ ^ 

increase up to 10 — 20 once rrih-^ decreases down to 40 GeV. Here we do not impose any 

experimental constraints, however, it is interesting that in fact only the points rrih^ > 90 

GeV can survive the LEP Higgs bound with the adopted choice of parameters. This is 

because the mixing with the doublet of order 10 % is necessary for mh 2 to reach 125 GeV. 

This region exactly matches with the area where the tuning is minimized. 

2 

For completeness, in figure 7 we plot A f and A^ 2 for Mq = 3 TeV and Mq = 5 

7?T 7YL Q 

TeV in addition to the Mq = 1.5 TeV case. The other parameters are fixed at tan/3 = 20 

777 ^ 

and = —100 GeV. We see that A^f is 10 — 20 for Mq = 3 TeV and ~ 40 for Mq = 5 

TeV. Also we see fine-tuning measure A |^2 ^ 10 if ruh^ > 60 GeV. It is remarkable that 
even the 5 TeV case is acceptable in the standard of the conventional models build around 
1 TeV. Again the LEP bound for hi is weak where the fine-tuning measures are minimized 
{nih^ ~ 100 GeV). 


5.3 Higgs couplings and collider signature 

In this section, we discuss the couplings of the Higgs bosons to the SM particles. They play 
crucial roles in production and decay of the Higgs bosons at the LHG and future colliders 
like ILG, FGG-ee(TLEP)/GEPG. The precision measurements of the signal strength for 
various production and decay channels will constrain the couplings of the SM-like Higgs 
boson at 0(1)% level or below [73-78]. The pattern of deviations from the SM predictions 
will bring hints for underlying models beyond the SM. Also there are chances to detect 
extra light Higgs bosons if their couplings to the SM particles are strong enough. 

The interaction of the GP-even Higgs bosons h, (i = 1 — 3) with the SM fields can be 
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H 

hs 


cos^ 2/3 -b ^ sin^ 2/3 

A\ + (m| — sin^ 2/3 

X^v‘^ sin^ 2/3 

/-il 

1 

0 

^ sin 2/3 

(-ii 

1 

— tan/3 

ir tan/3 

/-ii 

1 

cot f5 

o 


Table 1. Coupling constants of the CP-even Higgs bosons in the k = m| = 0 limit. 


described by the following effective Lagrangian [79]: 


£ = E 


2=1 


^ V ^ yj2v ^ •' yj2v 


cv cv 


TTV 


'V2' 


TTV 




(5.5) 


where we assume the custodial symmetry relating the W and Z couplings. The sum of / 
runs all the SM fermions. and are the field strengths of the gluon and photon, 
respectively. The coupling constants in the first line are normalized so that Cy^ = = 

1 for the SM Higgs boson hsM at the tree level. In the NMSSM, they are given by 




Sl^ sin/3-|-5jj cos/3, 


Iz 


■'2i 


cos /3 
sin/3 


if = e,n,T, d,s,b) 
if = u,c, t) 


(5.6a) 


(5.6b) 


where Sij is the diagonalization matrix for the mass matrix of the CP-even Higgs bosons. 
Cy is often quoted as in the literature and satisfies ~ 1 - have already 

depicted in the lower right panels in figure 1 and 2. These couplings are summarized 
in Table 1 in the k = = 0 limit. The singlet-like Higgs hs seems to be the lightest. 

However, once we include corrections from non-vanishing k and mg, it could be heavier 
than h for small tan/3. The light doublet h behaves like the SM Higgs boson within the 
approximation. Note that Cy and C!^ are stable against the small mixing between h and H 
induced by the finite n and m|, however, could be affected considerably by the tan/3 
enhanced [44, 49, 80]. 

The coupling constants C* and Cl^ are loop induced. Generally, if we have a scalar S, 
a fermion / and a (color-less) vector boson fields with the following couplings with the 
Higgs bosons [79], 




2=1 


-c: 


v/ 2 : 
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hiS^s - cj^hiff + 

■' y/SV ^ V 


^pip^ 


(5.7) 
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their contributions to Cg^ and (7^* are given by 


+ 2T{rf)C}Af{T}), 


N{rs)Ql i . (i. 

= - - - C,A,(t,) +-^- CfAfiTj) 


7Ql 


6 


CMri), 


(5.8a) 

(5.8b) 


where T{rs) (T(rf)) denotes the Dynkin index of the scalar (fermion) field in the repre¬ 
sentation Vs {vf) as defined by tr(T“(r)T^(r)) = T{r)5°‘^ with the SU{3) generator T°‘{r). 
N{rs) {N{rf)) is the number of freedom in the representation Vg {ff). Qs,f,p ^-re charges of 
the scalar, fermion and vector fields, respectively. The mass functions Asj^v{t) are defined 
as 


As{t) = ;^[/(t) -r], 

^/('t) = ^[(r - l)/(r) r], 

“ l)/('^) + 3r 2t\ 


f 

arcsin^ y/r 
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(5.9a) 

(5.9b) 

(5.9c) 


(5.9d) 


with r* = Since Agj^^ (r) decouple quickly for r >> 1, the standard model contribu- 

J m- ’ > 

tion is dominated by those of the top quark and the W boson and given by [79], 


For r << 1, we obtain 


= ^Afirf^) ^ 1.03, (5.10a) 

= lAfir,^^) - Ia,{t^^) « -0.81. (5.10b) 

Agj^v — 1- Then, we have. 


SClc:.'^Cl + 2T{rf)C}, 

N{rj)Q) 7Ql 

^ Cf - —Cg. 


(5.11a) 

(5.11b) 


Together with the relation. 
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(5.12) 
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we arrive at 
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^ 8\/2 


T{rs) 91ogdet(A^s) 4r(rj) 91ogdet(A4p 
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AN{rf)Q‘j Slogdet(Atj 
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(5.13a) 


(5.13b) 


where f ,, is the mass matrix for the corresponding held. Noticing that the coefficients of 
the derivatives are given by the steps of the beta function coefficients after integrating out 
the corresponding helds, this can be cast into the well known low energy theorem [81-83]: 


i as ^IdlogUsin) i a ^lSloga(/i) 

^12v/27rr; 4 dhi ’ V27rv 4 dhi 


(5.14) 


if all the integrated helds are much heavier than the renormalization scale, /i. Here, hi 
dependence of as and a enters through the mass thresholds. The contribution from the 
light helds should be included separately. 

In SUSY models, additional contribution comes from the charged Higgs boson, sfermion 
and chargino loops. The charged Higgs contribution to 5C^'^ is calculated as 




+ gl)s\.\i2l3 + \{gl+gl)cos^ 213} {i = h) 
48 2 

1 1 

< —{{-X^+gl)cos2l3 + -{gl+gl)W2P]^^ {i = H) 
48 2 mi 

— (2A + Ksin2,0)-^^ [i = S). 

24 mf _ 


(5.15) 


They are suppressed by gv/m'^_ and negligible in our model if m/j-(~ Mq) is 

above the TeV scale. 

The sfermion and chargino mass matrices are given by 


Ml 



( "ig, + vIaI + {vl - vj) f) yuiiAuiVu - yvd) 
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(5.16a) 

(5.16b) 


(5.16c) 


(5.17) 
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where Vu4 = {Hu,d)- The determinants of these mass matrices are calculated as 


det {Ml) = ~ ^ 


5? , gl 
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(-^5) = n + ydA {ml. + - {Adi “ /itan/3)^| 

2 ^ 

+ {vl - vl) { (f| + f) 
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2 ^ 

5'l+fi'2^ , ffl 2 1 , 4n 

j ""E? + J ^ ^ 


+ (vl - vj) 

det {M^-) = M 2 H - glvuVd- 


(5.18a) 


(5.18b) 


(5.18c) 

(5.18d) 


Thus the squark-slepton contributions to Ct and C\ are suppressed by lvn}_- at leading 
order and almost negligible if the squarks and sleptons are above the TeV scale. In addition, 
they are further suppressed if the following conditions are satished by the helds coupling 
via non-vanishing Yukawa interactions, 


“ i^ui - g cot /3)^ Ri 0, (5.19a) 

mg. + - {Ad^ - /xtan/3)^ ss 0, (5.19b) 

ml.+m% — (A 4 . —/xtan/3)^ Ri 0. (5.19c) 

It is interesting to observe that this is exactly the prediction of the TeV scale mirage 
mediation. 


CHu = 0 , = 1 , 

CHa = 1, = 1, Cf^ + Cf^ = 1, (5.20) 

with /r ~ Mq /tan /3 If any of the squarks or sleptons are light (e.g. = 0), this provides 

an exception of the anti-correlation between the hne-tuning of the EW symmetry breaking 
and the new physics contribution to the Higgs couplings [84]. 

The chargino contribution can be comparable to the SM contribution if the chargino 
mass is light enough [44, 79, 85, 86]. In the k = m| = 0 limit, it is calculated as 
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®This suppression also occurs in the type I model of [24]. 
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Table 2. Coupling constants of the CP-odd Higgs bosons in the limit k = 0, where tan 7 = Xv/Ax- 


For i = H we multiplied a suppression factor, As{t) ~ —12/ tan^/3. If /r 'D ^ \ r-^ 1 and 
once the finite k and induce a mixing of 0 ( 10 )% between the singlet and the heavy 
doublet, then 0 ( 10 )% deviation is expected in O^. 

The effective Lagrangian describing the CP-odd Higgs bosons is given by 
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(5.22) 


The renormalizable couplings between the CP-odd Higgs boson and the vector bosons are 
forbidden by the CP symmetry. We omit the loop induced couplings with W/Z bosons 
because they do not give sizable effects in the production and decay of the Higgs bosons. 
The coupling constants with the fermions are given by 


c- = <: 


(pt) 

(pt) 


cos /3 
sin (3 


if = e,n,T,d,s,b) 

if = u,c, t) 


(5.23) 


where V is the diagonalization matrix for the mass matrix of the CP-odd Higgs bosons. 
These coupling constants in the k = 0 limit are summarized in the table 2. Note that those 
of ai are highly suppressed by sin 7 ~ Xv/A\ because they require the mixing with the 
heavy doublet. The coupling constants with gluon and photon are calculated as [87, 88 ], 

K (5.24a) 

= (5.24b) 

where the mass function is given by 

4(t) = 2T-V(r). (5.25) 

In the limit << 1, they are approximated as 

5Cl T{rf)C'/, (5.26a) 

6C\ ~ N{rf)Q)C'f\ (5.26b) 
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Figure 8. The scale factors of the CP-even Higgs coupling constants for tan/3 = 3, Mq = 1500 
GeV and = —100 GeV. The vertical lines show the expected precisions of the SM-like Higgs 
coupling constants for LHC and ILC, A(kx)^ « 2Akx- 


The charged Higgs contribution is forbidden by the CP symmetry and the sfermion con¬ 
tribution is highly suppressed as in the CP-even case. The chargino contribution to is 
approximated as 


7 y 


a=l 


- = 2 
Xa 




iJ- 


4-i3 


sin /3 cos^ j3 


\v 
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^^(lntaiH/3-/7r)^ 


(/ = 1 ) 
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(5.27) 


For i = 2, we multiplied a suppression factor H((t)/2 inside the bracket. If /r ~ n, A ~ 1, 
could be comparable with . 

Next we present the results of our numerical calculation on the coupling constants in¬ 
cluding the effect of finite k and rrig, the radiative corrections and various phenomenological 
constraints. Here we use the scale factors of the Higgs couplings relative to the SM instead 
of the coupling constants themselves [89, 90], 


'‘X = (5.28) 

In figure 8, we plot the k’s for tan/3 = 3, Mq = 1500 GeV and = —100 GeV. 
The horizontal axis is chosen as the mass of the singlet-like Higgs boson which is heavier 
than the SM-like Higgs boson for small tan/3. We impose the condition 124GeV < ruh^ < 
126 GeV and various phenomenological constraints built in the NMSSM Tools package. 
The parameter region with 124 GeV < rrih^ <126 GeV is excluded by the LEP Higgs boson 
search [36]. The vertical lines show the expected precision of the coupling constants for 
LHG (300/6-1) [ 73 ]^ jlq (250 GeV) and ILG (250 GeV+500 GeV) [74] (The position of 
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Facility 

LHC 

HL-LHC 

ILC(250) 

ILC(500) 

ILC(IOOO) 

ILC(LumUP) 

xA (GeV) 

1,400 

1,400 

250 

250+500 

250+500+1000 

250+500+1000 

J£dt (/b-i) 

300/expt 

3000/expt 

250 

250+500 

250+500+1000 

1150+1600+2500 


5-7% 

2-5% 

17% 

8.3% 

3.8% 

2.3% 

Kg 

6-8% 

3-5% 

6.1% 

2.0% 

1.1% 

0.7% 

Kw 

4-6% 

2-5% 

4.7% 

0.4% 

0.3% 

0.2% 

KZ 

4-6% 

2-4% 

0.7% 

0.5% 

0.5% 

0.3% 

He 

6-8% 

2-5% 

5.2% 

1.9% 

1.3% 

0.7% 

Hb 

10-13% 

4-7% 

4.7% 

1.0% 

0.6% 

0.4% 


14-15% 

7-10% 

6.4% 

2.5% 

1.3% 

0.9% 


Table 3. Expected precisions on the Higgs coupling for LHC [73] and ILC [74]. Two numbers 
for LHC represent different assumptions on the systematic error. The seven parameter HXSWG 
benchmark parametrization [89, 90] is assumed. 


the center is arbitrary). We also quote the numbers in table 3. We plot the precision of k\ 
for Ky and take conservative values for the LHC expectation. 

The pseudo Nambu-Goldstone nature of the SM-like Higgs boson well suppresses the 
mixing with the singlet and the deviation from the SM is at most 20 %. This suppression 
is also important for not reducing the lighter SM-like Higgs boson mass due to the mixing. 
The approximate sum rule = 1 means that the deviation mainly comes from 

the mixing with the singlet and the SUSY contribution is negligible except for h (r) for 
which the contamination of the heavy doublet has tan/3 enhanced coupling. The singlet¬ 
like Higgs boson is hidden from the collider search and an interesting target for the precision 
measurement at the LHC and ILC. In this parameter region, the prediction does not change 
significantly even if we adopt a conservative bound, 120 GeV < miy < 130 GeV, taking into 
account the theoretical error [57]. 

In figure 9, we plot the same scale factors (normalized by the SM Higgs couplings) for 
tan/3 = 10. Increasing tan/3, the singlet-like Higgs boson is now lighter than the SM-like 
Higgs boson and we impose 124 GeV < < 126 GeV. We again choose the horizontal 

axis as the singlet-like Higgs mass. The tree-level contribution to the SM-like Higgs mass 
is now not effective with larger tan /3 and small mixing with the singlet is favorable to 
achieve 125 GeV. However, a strong constraint exists for the light singlet from the LEP 
Higgs boson search if it couples with the SM particles through the mixing with the SM-like 
Higgs boson [36]. The khaki shaded region is mostly excluded by this constraint. The 
region 90 GeV < < 100 GeV survives due to a small excess in the LEP measurement 

[36, 40, 54-56]. Around rrih-^ ~ 120 GeV, the two bosons almost degenerate and the mixing 
is maximized, but the LEP bound is weak in this region. For large A, the new radiative 
corrections to the SM-like Higgs mass, eq.(4.11) could relax the mixing with the singlet 
(around ruh^ ~ 40 GeV). In this parameter region, the lightest neutralino becomes lighter 
than the half of the Z boson mass and the decay channel, Z —)• 2xi opens. Then the region 
is excluded by the measurement of the Z boson width [91]. The mixing with the singlet 
leads to large deviations in the coupling constants of the SM-like Higgs boson observable at 
the LHC. On the other hand, searching the light singlet boson missed by LEP is important 
physics at ILC. The sum rule holds again except for 6(r) and g, for which the mixing with 
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Figure 9. The scale factors of the CP-even Higgs coupling constants for tan/3 = 10, Mq = 1500 
GeV and = —100 GeV. The vertical lines show the expected precisions of the SM-like Higgs 
coupling constants for LHC and ILC, A(kx)^ ~ 2Akx- and ( k ^)^ ((/^ v )^ ('^t )^) ™ the 

upper (lower) left panel almost overlap each other. 


the heavy doublet is not negligible due to the tan /3 enhancement. The small excess in 
around m/n — 30 GeV is contribution from the chargino interfering with that of the SM. The 
surviving region expands once we raise the SUSY breaking scale Mq and increase the top 
radiative correction to the SM-like Higgs mass. Also note that the prediction is sensitive 
to the theoretical error in the Higgs mass calculation. If we take a conservative bound, 
e.g. 120 GeV < < 130 GeV, much lighter hi (m/^^ > 20 GeV) is allowed because the 

singlet-doublet mixing is suppressed. For such the hi couplings to the SM particles 

and the deviations of h 2 couplings from those of the SM Higgs boson are highly suppressed 
except for b and r. 

In hgure 10, we show the scale factors of the lightest GP-odd Higgs boson for the 
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Figure 10. The scale factors for the CP-odd Higgs boson coupling constants for tan /3 = 3 (upper 
panel), tan/3 = 10 (lower panel), Mq = 1500 GeV and = —100 GeV. 


same parameter set with the CP-even bosons. The horizontal axis is chosen as the mass 
of ai, which is rather sensitive to the choice of the PQ symmetry breaking parameter, A^- 
and (k^^)^ are doubly suppressed by tan“^/3 and sin^y and negligible. The direct 
production of oi at the LHC is almost hopeless and the gamma-gamma option of ILC may 
help if (kIj,^)^ is large enough. 

5.4 Dark Matter 

The existence of the dark matter (DM) is one of the few compelling evidences of physics 
beyond the SM. In SUSY, the Lightest Supersymmetric Particle (LSP) is stable and good 
candidate of the dark matter, once we assume the R-parity to forbid the dimension 4 
baryon/lepton number violating operators to prevent rapid nucleon decay. If LSP is weakly 
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interacting massive particle (WIMP), its annihilation cross section predicts the thermal 
relic abundance in the ballpark of the observed DM density. This ’’WIMP miracle” scenario 
provides a valuable constraint on the SUSY models in light of the recent progress of precision 
cosmology (for more details, see e.g. [92]). However, in the mirage mediation, the oscillation 
of the modulus field after the inflation once dominates the energy density of the universe. 
The reheating temperature of the modulus field is 0(100) MeV. It is much larger than the 
temperature at which the Big Bang Nucleosynthesis (BBN) takes place, however, lower than 
the decoupling temperature of LSP. Huge entropy produced by the modulus decay dilutes 
the preexistent thermal LSP abundance to negligible level. The modulus decay produces 
non-thermal LSP. In particular, daughter gravitino later decays at the decay temperature 
0(10) MeV and overproduces LSP (moduli-induced gravitino problem [94, 95]). This non- 
thermal overproduction of DM does not mean the model is excluded, if we consider a specific 
non-minimal cosmological scenario or modify the DM sector of the model. For example, 

1. Suppressing the modulus oscillation or reducing the gravitino branching ratio (e.g. 
enhanced symmetry [96], adiabatic oscillation [97, 98], alignment [99]). 

2. Violating R-parity and introducing another source of DM (e.g. axion [100]) or intro¬ 
ducing thermally decoupled light LSP (e.g. axino [31]). 

3. Diluting the LSP after the gravitino decay (e.g. by thermal inflation [101-103], Q-ball 
[104] or unstable domain wall [105]). 

Explicit model building is beyond the scope of this paper. In this section, we just present 
the prediction of our model with the minimal DM sector assuming the first or the third 
scenario. The calculation is performed based on micrOMEGAs built in NMSSMTools [64]. 

In our model, LSP is given by mixture of the higgsino and singlino. Generally, the 
thermal relic of the pure higgsino of 0(100) GeV cannot saturate the observed DM density. 
However, mixing with the singlino reduces the annihilation cross section and enables to 
explain the observed amount of DM (Well-tempered neutralino) in the first scenario as in 
the case of bino and higgsino in MSSM [40, 106]. (There is no constraint in the third 
scenario.) 

The searches of DM are classified into the direct and indirect detections. The former 
detects the scattering of target nucleus by DM at underground laboratories. The latter 
observes the annihilation/decay products of DM in celestial/astrophysical circumstances 
such as the center of the earth/sun or galaxies. 

We first discuss the prediction for the direct detection. In the direct detection ex¬ 
periments, the scattering between (slow) neutralino and target nucleus proceeds through 
the scalar or pseudo-vector interactions due to its Majorana nature. The former gives the 
spin-independent cross section and the latter gives the spin-dependent cross section [92]. 
In our model, the scalar interaction mainly comes from the exchange of the GP-even Higgs 
bosons and the pseudo-vector interaction mainly comes from the exchange of the Z boson. 
The spin-independent nuclear cross section is enhanced by the square of the mass number 
because the process is coherent, while the spin-dependent cross section does not have such 
a factor. Thus the experiments are much sensitive to the spin-independent cross section 
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Figure 11. The predictions for the direct and indirect searches of the dark matter for tan/3 = 3, 
Mo = 1500 GeV and = -100 GeV. 


than the spin-dependent ones at nucleon level. On the other hand, the Z boson interacts 
with nucleon much stronger than the Higgs bosons. Therefore both cross sections could 
give meaningful constraints on our model, depending on the amount of higgsino component 
in LSP. 

In the upper left panel of figure 11, we show the spin-independent cross section of LSP 
for tan/3 = 3 case. The red region indicates the region where the lightest CP-even boson 
corresponds to the SM Higgs boson and the blue region indicates the region where the 
second lightest CP-even boson is the SM-like Higgs boson (although excluded by the other 
phenomenological constraints). In the pink region, the thermal abundance saturates the 
observed value of DM density. Around misp ~ 100,350 GeV, the scattering amplitudes 
mediated by the lightest and second lightest CP-even bosons cancel each other and the cross 
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section vanishes [37, 107]. The neutralino-Higgs boson coupling involves A, therefore the 
red region is on the upper edge of the distribution because almost maximum A is required 
to raise the Higgs mass to 125 GeV. With the particular model parameters chosen in the 
hgure, the red region overlaps with the pink region where the thermal abundance saturates 
the DM, however, is almost excluded by the LUX experiment [108]. If we increases Mq, 
the cross section decreases while the pink region remains almost same. In the upper right 
and lower left panels, we present the spin-dependent cross sections for proton and neutron, 
respectively. The red region is again on the edge of the distribution because the singlino- 
higgsino mixing terms in the mass matrix are proportional to A. The current experimental 
bounds are 1-2 order above the model prediction [109, 110]. 

Next we discuss the indirect detection. We consider the constraint by the gamma ray 
from the dwarf spheroidal galaxies observed by Fermi satellite [111]. The annihilation of the 
neutralinos mainly produces W~^W~, pairs, tt, bb pairs, hittj and hiZj pairs if they are 

kinematically available. Other channels are strongly suppressed due to the Majorana nature 
of the neutralino and CP conservation [92]. In the lower right panel of hgure 11, we show 
the total annihilation cross section of the neutralino and the experimental bound assuming 
the decay proceeds through W~^W~ and bb pairs. In our model, ti decay dominates the 
process above the threshold. The bound is one order larger than the prediction assuming 
that the gamma ray spectrum from ti is not much different from those of WW, bb pairs. 

Here we impose 124 GeV < rrih-^ < 126 GeV, however, the conservative Higgs mass 
bound, 120GeV < mhi < 130GeV does not change the results for tan/3 = 3 signihcantly. 

In hgure 12, we present the same set of the predictions for tan /3 = 10 case. The blue and 
cyan regions indicate the region where the second lightest CP-even boson corresponds to the 
SM Higgs boson. The blue region is already excluded by the phenomenological constraints 
explained in the previous section. In the pink region, the thermal LSP abundance saturates 
the observed DM density. In the surviving cyan region, the annihilation cross section is 
larger than the value saturating the observed DM density. The spin-independent cross 
section (upper left panel) shows a dip around misp ~ 100 GeV due to the cancellation 
between the amplitudes mediated by the two CP even bosons. While the spin-dependent 
cross sections (upper right and lower left panels) do not show such a behavior because only 
one diagram dominates the process. They are enhanced with tan /3 because fi ~ Mq /tan /3 
and the higgsino component in LSP increases with tan j3. The cyan region is marginally 
excluded by the spin-independent cross section and the spin-dependent cross section for 
neutron. Discounting the constraints considering the astrophysical, hadronic and nuclear 
physics ambiguities in the calculation (see e.g. recent analysis in [93]), it is interesting that 
the surviving region is exactly located around the dip structure where the hne-tuning of 
the EW symmetry breaking is also minimized. Once observed, this feature can be checked 
by the experiment using the target nucleus insensitive to the spin-dependent cross section 
(enriched even-even nucleus) or combining the results of different target nuclei. In the lower 
right panel, the resonance structure in the annihilation cross section appears due to the 
pseudo-scalar mediated bb decay, however, this region is excluded by the phenomenological 
constraints. Above the W~^W~, Z^Z^ threshold, the process is dominated by W^W~ and 
Z^Z^ decays. The surviving region is about a factor 2 below the current bound. Once we 
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Figure 12. The predictions for the direct and indirect searches of the dark matter for tan/3 = 
lO.aAAMo = 1500 GeV and = -100 GeV. 


raise the SUSY scale, Mg, the constraint from the LEP Higgs boson search becomes weaker 
until the radiative correction overshoots 125 GeV and the cyan regions in the hgures expand, 
while the EW hne-tuning deteriorates as shown in hgure 7. The conservative Higgs mass 
bound, 120 GeV < mh^ < 130GeV opens a lighter LSP region {tulsp — 50 GeV). However, 
this region is excluded by the upper bound of the spin-independent cross section, if the 
dominant component of the DM consists of the lightest neutralino. 

In this section, we discussed the DM prediction of the model assuming the LSP saturates 
the observed DM abundance. We stress that the novel structure of the EW symmetry 
breaking in our model is not directly connected to the DM physics and the model will not 
be excluded even if the prediction of the minimal model is not conhrmed by the experiments. 
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6 Conclusion 


We have studied the NMSSM with the TeV scale mirage mediation. Our choice of the 
modular weights (4.1) realizes a little hierarchy — Mq/Stt^ << Mq at the EW scale 
and automatic cancellation between and m|^^/tan^/3 in (1.1). This leads to a natural 
EW symmetry breaking with 0(10)% tuning for the TeV scale stop mass. The /x term 
can be as large as y/mzMo without introducing a hne-tuning despite conventional wisdom 
tells that it lies around the EW scale. The Higgs boson mass is lifted to 125 GeV by the 
new quartic coupling (small tan/3) or small mixing with the singlet (moderate tan/3) in 
addition to the radiative corrections (from the top and singlino sectors). The suppression 
of the singlet-doublet mixing is important for raising the Higgs boson mass in the hrst case 
and hiding the singlet from the LEP Higgs boson search in the second case. In the NMSSM, 
if K = 0 and = 0, the doublet Higgs VEVs spontaneously break an approximate scale 
symmetry (4.12). The resultant pseudo Nambu-Goldstone boson (mass eigenstate) almost 
consists of the doublets and corresponds to the SM-like Higgs with a suppressed singlet 
component. There is another scale symmetry if k = 0 and = 0 in which S and are 
interchanged. The corresponding Nambu-Goldstone boson is the singlet-like Higgs boson 
with suppressed component. The mixing with the suppressed components must break 
the two symmetries and disappear if Hd decouples. Thus it is highly suppressed although 
the symmetries are broken at the EW scale. To evaluate the hne-tuning, we introduced 
a simple formula for the hne-tuning measure having a clear physical meaning which calls 
for careful treatment of the radiative corrections. Using it, we performed a comprehensive 
analysis of the hne-tuning in the EW symmetry breaking which shows the model is tuned 
at most 0(10)% for Mq = 1.5 TeV. For moderate tan/3, the least tuned region is realized 
with the singlet Higgs boson mass around 100 GeV which overlaps with the anomaly found 
in the LEP Higgs boson search. 

We further investigated the coupling constants of the Higgs bosons in our model. The 
radiative corrections from the SUSY particles and heavy Higgs bosons are highly suppressed 
due to the little hierarchy. Thus the deviation from the SM occurs through the suppressed 
mixing between the SM-like and singlet Higgs bosons. The numerical calculation showed 
that 20 % deviations are possible in the squares of the scale factors, {nx)‘^ 
rule (k^)^ -|- (k^)^ = 1 holds well except for those of b and r. The mixing with the heavy 
doublet is also important for b and r due to tan/3 enhancement and an 0(1) deviation is 
possible for their couplings to the SM-like Higgs boson. Future precision measurements at 
the LHG and the ILG will provide good opportunities to explore these deviations and the 
hidden singlet-like scalar. 

We also studied prospects for the direct/indirect dark matter searches assuming the 
minimal dark matter sector in our model. The spin-independent neutralino-nucleus scatter¬ 
ing is mainly mediated by the SM-like and singlet-like Higgs bosons. The two amplitudes 
interfere destructively in some parameter regions. This suppresses the cross section below 
the current experimental bound. While the spin-dependent scattering is dominated by the 
Z boson exchange and does not show such a feature. The surviving regions from various 
phenomenological constraints are almost excluded for tan /3 = 3 and on the verge of exclu- 
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sion for tan /3 = 10. The suppression of the spin-independent cross section in the surviving 
region is characteristic signature to identify the model. The bound on the neutralino an¬ 
nihilation cross section from the gamma ray emission by the dwarf spheroidal galaxy is an 
order above the prediction and will soon start to explore our model. 


A Soft SUSY breaking terms 


Here we give explicitly soft SUSY breaking terms induced by the mirage mediation mech¬ 
anism in the NMSSM. 

In the mirage mediation, the soft parameters at the scale just below Mqjjx are given 

by 


MaiMcur) = Mo + 

^ 3/2 

Aijk{MQUT) = (Q + Cj -|- Ck)Mo — (ji -|- Jj -|- 7A:) ’ 

mfiMcur) = (^- 1 ) 


where 


i 

a jk 

a jk 


ii = 


0,2 d-fi 
din HR 


(A.2) 


Here, T^(Adj) and T^((/)*) denote Dynkin indices of the adjoint representation and the 
representation of matter fields 0*. We have assumed Uij = ’^kiyijkVjki to be diagonal. 

Within the framework of the NMSSM, the /3-function coefficients, anomalous dimen- 
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sions and other coefficients in the above equations are obtained as 


bs = -3, 62 = 1, 61 = 11, 
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So 1 9 ^ 9 

= -52 + 2^1 - ^ ’ 

75 = -2t^ - 2\\ 

iQa = + Ygfl'l “ “*■ yb)b3a, 

3 2 3 2 O 2 r 

lUa = ^53 + g 5 l “ ^Vt bSa, 

82 ^2 n 2 r 

7 d, = ^53 + g 5 l - 25 fe() 3 a, 

7La = ^gi + ^5? - y%a, 

= 2gl - 2y%a, (A.3) 


= 3^2 + 5i - ^ytC'HuQiUI - 

OHa = 352 + 5 l - ^yl^HiQzDl “ 2 y^a/f^I, 3 E| - 2 X^aH^HaS, 

6s = -2\^aH^HdS - I^^asss, 

&Qa = Y^l + 352 + g5i - 2(5t a/f„Q3t/| + ylaHaQiDi)b3a, 

16 2 16 2 3 2 e 

(^Ua = -1793 + -^9l - maH„Q3f/|<^3a, 

16 2 4 2 . 2 

= Y ^3 + g 5 l - ^yb(^HaQ 3 D^ 03 a, 

Ola = 3^2 + 5l - ‘^ylaH^LsE^ha, 

Osa = ^gl - (A.4) 


= \92 + y5f - ^ytbyt - A^6 a, 
iHa = ^gi + y9i ~ ^ybK - y%T - a^^a, 

75 = - 2 k^ 6 k - 2 A^ 6 a, 

3 11 

7Qa = “3^3 + -5(| + —gj — {utby^ + 5fe6yJ(J3a, 

88 

iUa = -853 + y5f - “^VibytSsa, 

22 

7D„ = -85I + y5f - ‘^Vbby.bsa, 

3 4 11 4 27 C 

7L. = 25^2 + yS'l - Vr^yMa. 

7s„ = 22fft - 22/>y,<53a, (A.5) 
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where 


K = -yffi - 3^1 - -jg! + 6y? + vl + A^, 

Kb = “y^S “ 3^2 - g5l + y* + %6 + 2/r + A^, 

= —35^2 — ^gi + 3yb + 4y^ + A^, 

6k = 6A^ + 6k^, 

6a = —35^2 ~ 51 + 35t + 35b + 5 t + ^A^ + 2k^. (A.6) 

Here, Qa, Ua, Da, La, and Ea denote left-handed quark, right-handed up-sector quark, right- 
handed down-sector quark, left-handed lepton, and right-handed lepton helds, respectively, 
and the index a denotes the generation index. We have included effects due to Yukawa 
couplings, yt,yb, and y^, only for the third generations. 
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